Presupernova collapse models with improved weak-interaction rates 
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Improved values for stellar weak interaction rates have been recently calculated based upon a large 
shell model diagonalization. Using these new rates (for both beta decay and electron capture), we 
have examined the presupernova evolution of massive stars in the range 15-40 M©. Comparing our 
new models with a standard set of presupernova models by Woosley and Weaver, we find significantly 
larger values for the electron-to-baryon ratio at the onset of collapse and iron core masses reduced 
by approximately O.IMq. The inclusion of beta-decay accounts for roughly half of the revisions, 
while the other half is a consequence of the improved nuclear physics. These changes will have 
important consequences for nucleosynthesis and the supernova explosion mechanism. 



PACS numbers; 26.50.+X, 97.60.Bw, 97.10. Cv 

The late stages of massive stellar evolution are strongly 
influenced by weak interactions which act to determine 
the core entropy and electron-to-baryon ratio, Yg, of the 
presupernova star. Electron capture reduces the number 
of electrons available for pressure support, while beta de- 
cay acts in the opposite direction. Both processes gener- 
ate neutrinos which, for densities p < 10^^ g/cm'^, escape 
the star carrying away energy and entropy from the core. 

Electron capture and beta decay during the final evo- 
lution of a massive star are dominated by Fermi and 
Gamow- Teller (GT) transitions. While the treatment 
of Fermi transitions (important only in beta decays) is 
straightforward, a correct description of the GT transi- 
tions is a difficult problem in nuclear structure. In their 
pioneering work on the subject Fuller, Fowler and New- 
man (FFN) Q estimated electron-capture rates assuming 
a single GT resonance. The properties of this resonance 
were derived on the basis of the independent particle 
model, supplemented by Fermi contributions and experi- 
mental data for low-lying transitions, whenever available. 
These authors also noted the importance of the 'backres- 
onances' for beta decay. These are excited states in the 
decaying nucleus which are connected by strong GT tran- 
sitions to low-lying states in the daughter nucleus and, by 
thermal population and with increased phase space, can 
significantly contribute to the stellar beta decay rates. 

Recent experimental data shows that the GT distri- 
butions in nuclei are quenched, compared to the inde- 
pendent particle model value, and strongly fragmented 
over many states in the daughter nucleus. Both effects 
are caused by residual interaction among the valence nu- 
cleons and an accurate description of these correlations 
is essential for a reliable evaluation of the stellar weak 
interaction rates due to the strong phase space energy 
dependence, particularly of the stellar electron capture 
rates. The shell model is the only known tool to reli- 
ably describe GT distributions in nuclei ||^. Its appli- 
cation to iron mass nuclei in the middle of the p/-shell 
as required in the presupernova collapse, however, has 
long been inhibited due to the extremely large model 



space dimensions involved. After significant progress in 
shell-model programming |^ and hardware development 
the situation has changed very recently and in Ref. Q| 
it has been demonstrated that state-of-the-art diagonal- 
ization studies, typically involving a few 10 million con- 
figurations, are indeed able to reproduce all relevant in- 
gredients (GT± strength distributions for changing pro- 
tons (neutrons) into neutrons (protons), level spectra and 
half-lives) and hence have the predictive power to reliably 
calculate stellar weak interaction rates. This program has 
recently been finished and stellar weak-interaction rates 
for nuclei with A — 45-65 have been calculated based 
on the shell-model results, supplemented by experimen- 
tal data, wherever available. The shell-model rates have 
been discussed and validated in It has been found 
that for p/-shell nuclei the shell-model electron-capture 
rates are smaller than the FFN rates by, on average, an 
order of magnitude, for the reasons explained in ||^. The 
situation is different for the beta decay as the shell model 
and FFN rates are of the same magnitude for the most 
relevant nuclei to be identified below. 

To study the influence of the shell model rates on pre- 
supernova models we have repeated the calculations of 
Woosley and Weaver |^ (henceforth WW) keeping all 
the stellar physics, except for weak rates, as close to the 
original studies as possible. The present calculations have 
incorporated the new shell-model weak interaction rates 
(including electron capture, positron emission, and beta 
decay) for nuclei with mass numbers A = 45-65, supple- 
mented by rates from Oda et al. ^ for lighter nuclei. In 
practice, the weak rates for these lighter nuclei were not 
very important for determining the presupernova struc- 
ture, but only dominate prior to silicon burning. We 
also note that, for sd-shell nuclei with A — 17 — 39, the 
FFN rates agree rather well with the shell model rates 
previously determined by Oda et al. 0| . 

The earlier calculations of WW, to which we compare, 
used the FFN rates for electron capture and an older set 
of beta-decay rates taken from ||] and |^ . Shortly after 
the models of WW were calculated, it was recognized 
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that these older beta-decay rates were inadequate and 
that use of the larger values from FFN would appreciably 
alter the results These are the first models of the 

WW variety to incorporate realistic beta-decay rates and 
electron-capture rates in a complete stellar model, even 
though beta decays were included in some other models 
of massive stellar evolution 1 11 1 . In a separate paper [n2| , 



we will present the comparison of models that use the full 
FFN rate set and the new rate set, and in a future paper 
we will examine the evolution of stars of other mass and 
metallicity. 

Fig. |l| shows the late evolution, following core oxy- 
gen burning, of the central temperature and density in 
15 M0 and 25 stars as well as the central value of 
the electron-to-baryon ratio. Ye- Time is measured here 
backwards from the time of iron core collapse, which is ar- 
bitrarily t = 0. Silicon burning ignites with a mild "flash" 
and the core becomes convective and expands. For the 
25 M0 star the temperature and density trajectories in 
the new and old models are similar, but the calculations 
with shell model rates have significantly larger values of 
Yg. This difference persists throughout the iron core, not 
just at its center. Larger values of Ye also result in the 
15 Mq star, but there even the density and temperature 
structures of the presupernova star are appreciably al- 
tered. 

To understand the origin of these differences, we will 
now explore the role of the weak-interaction rates in 
greater detail. Weak processes become particularly im- 
portant in reducing Yg below 0.50 after oxygen depletion 
(~ 10^ s and 10^ s before core collapse for the 15 Mq and 
25 Mq stars, respectively) and Y^ begins a decline that 
becomes precipitous during silicon burning. Initially elec- 
tron capture occurs much more rapidly than beta decay. 
Since the shell model rates are generally smaller than 
the FFN electron capture rates, the initial reduction of 
Ye is smaller in the new models. The reduction is less 
pronounced during silicon shell burning both because the 
evolution time scale becomes quite short and because nu- 
clei near the valley of beta-stability, with smaller weak- 
interaction rates, have already been produced (from a 
composition that initially had N = Z). During this pe- 
riod the core matter is composed of nuclei with yl < 65 
which are carried in the calculation. 

An important feature of the new models is that beta 
decay becomes temporarily competitive with electron 
capture after silicon depletion in the core and during 
silicon shell burning. That this would occur was fore- 
seen by jl^ on the basis of one-zone models. Here we 
see it occurring in a complete stellar model. Moreover, 
the new electron capture rates are smaller than FFN and 
thus offer less resistance to beta decay. Dynamic weak 
equilibrium, in the sense described by | [l0[ |, thus occurs 
at larger values of Yg. Interestingly, by the time the iron 
core is actually collapsing, weak equilibrium no longer 
exists. The increase in density closes the phase space for 



beta decay and electron capture again predominates by 
a large factor. It is this special characteristic of "presu- 
pernova models" that lead some researchers in the past 
to miss the importance of beta decay during a transient 
stage an hour or so prior to collapse. 

While dynamic weak equilibrium is achieved in the 
15 M0 model, with the new rates, it is not in the 25 Mq, 
though beta-decay still offers a non-negligible resistance 
to electron capture even there. This is in part due to 
the shorter time scale of silicon shell burning in the more 
massive star and also the larger value of Ye in the cores 
of stars with higher entropies (Fig. |^ . 

The presence of an important beta-decay contribution 
has two effects. Obviously it counteracts the reduction of 
Ye in the core, but equally important, beta decays are an 
additional neutrino source and thus add to the cooling 
of the core and a reduction in entropy |lO|. This cool- 
ing can be quite efficient, as often the average neutrino 
energy in the involved beta decays is larger than for the 
competing electron captures. As a consequence the new 
models have significant lower core temperatures than the 
WW models after silicon burning, which is particularly 
pronounced for the 15 Mq star. During the contraction 
stage electron capture is again more important than beta 
decays, associated with the increased electron Fermi en- 
ergy. Although the shell model rates are individually 
smaller than the FFN electron capture rates, the effec- 
tive electron capture rate is larger in the new models 
as the evolution now proceeds along a trajectory with 
larger Ye values involving nuclei with smaller Q-values 
thus making electron capture energetically easier. 

In summary, the shell model weak interaction rates 
result in significantly larger Ye values during the presu- 
pernova evolution of a 15 and 25 Mq star than calculated 
by WW. Part - about half - of the change is due to in- 
cluding beta-decay and the other half is due to slower 
rates for electron capture. Fig. ^ shows that this is a 
general finding for all stars in the mass range 11-40 Mq. 
The central values of Ye at the onset of core collapse are 
increased by 0.01-0.015. This is a significant effect. We 
note that the new models also result in lower core en- 
tropies for stars with M < 20 Mq. The larger core Ye 
values and the lower entropies suggest that these stars 
will have a larger homologous core than currently as- 
sumed. For AI > 20 Mq, the new models actually have 
a slightly larger entropy. Thus one might expect that in- 
creased electron capture rate on (the more abundant) free 
protons will partly counteract the increase in Ye values 
during the subsequent core collapse phase. In general, 
core collapse calculations with detailed neutrino trans- 
port are required before definite conclusions about the 
explosion mechanism can be drawn. We will provide our 
models to those wishing to attempt such calculations. 

Another important property for the core collapse is 
the size of the iron core, which we define as the mass 
interior to the point where the composition becomes at 
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Time till collapse (s) Time till collapse (s) 

FIG. 1. Comparison of the time evolution of key quantities at the center of a ISM© (left) and 25 M© (right) star between 
the models by Woosley & Weaver models (WW, dashed) and the present ones using the shell model weak interaction rates 
(LMP, full lines). The upper panels show the temperature (thick lines) and density (thin lines), the middle panels show the 
electron-to-baryon ratio, Y^, and the lower panels show the time-derivative of Ye due to electron capture and beta decay for 
the models using the LMP rates. As reference points, in the ISM© (25 M©) star central silicon burning ignites at 7 x 10^ s 
(5.6 X 10* s), terminates at 6.6 x 10* s (10* s), silicon shell burning starts at 3.2 x 10* s (4.7 x 10^ s) and core contraction set in 
at t 280 s (~ 230 s) before core collapse. Note that the dYf./dt shown in the lower panels includes only the contribution of 
weak interaction processes to the total change of the central Ye; during central silicon burning convection mixes down matter 
with higher Ye from layers above |ia|. 



least 50% of iron group elements {A > 48). As is shown 
in Fig. H, the iron core masses are generally smaller in 
the new models. The effect is larger for stars more mas- 
sive than 20 Mq , while for the most common supernovae 
{M < 20 Mq) the reduction is about 0.05 Mq. This re- 
duction appears to be counterintuitive at first glance with 
respect to the slower electron capture rates in the new 
models. It is, however, related to changes in the en- 
tropy profile during silicon shell burning which reduces 
the growth of the iron core just prior to collapse [p^ . 
Clearly, though we have, for reasons of space concen- 
trated upon central values, it is the entire distribution of 
entropy and Ye in the stellar interior that determines its 
final evolution. 

Confidence in the shell model rates stems from the 
recent measurements of GT distributions in iron mass 
nuclei which are all well reproduced by the shell model 
calculations |Q. However, the energy resolutions of these 
pioneering (n,p) charge-exchange studies performed at 
TRIUMF has been rather fimited (~ 1-1.5 MeV) and 
they have been performed for stable nuclei only. These 
limitations are likely to be overcome in the near fu- 
ture as measurements with charge-exchange reactions like 
(d, ^He), (t, ^He) promise data with an order of magni- 
tude improved resolution. Furthermore, after radioac- 



tive ion beam facilities will become operational it will 
be possible by inverse techniques also to determine the 
GT distributions for unstable nuclei. Of course, lab- 
oratory experiments cannot measure directly the rele- 
vant stellar rates as these involve, for example, electron 
capture or beta decays from excited states. Neverthe- 
less, high-resolution charge-exchange or beta-decay ex- 
periments are important for two reasons: First, they are 
stringent constraints for the nuclear models and their pre- 
dictive powers. Second, such experiments can determine 
the energy positions of the daughter states for the GT 
transitions which can then be used directly in the deter- 
mination of the rates. 

To guide such experiments we have attempted to iden- 
tify the most important nuclei for electron-capture and 
beta decay during the final stages of stellar evolution. 
The relevant quantity is the product of abundance of the 
nuclear species in the core composition and electron cap- 
ture rate (or beta decay rate). Table || lists the most 
important nuclei at selected points during the final evo- 
lution of our 15 Mq and 25 Mq stars. Since the Ye values 
are significantly larger in the new models, the important 
fiows now involve nuclei much closer to stability so that 
several of the most important electron capturing nuclei 
(e.g., 54,56,58pg^ ^^Mn, ^^Cr) are stable. Beta decay, on 
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TABLE I. Most important nuclei for electron capture and beta decay at selected points (characterized by temperature T, 
density p and electron-to-baryon ratio Ye) during the final evolution of a 15 M0 and 25 M0 stars. The total electron capture 
Acc and beta decay A^- rates are listed as well as the 3 dominating nuclei; the number in parentheses defines their percentage 

to the respective rates. 



T(K) 


P (g cm ^) 




Aec (s-1) 




electron capture 


A^- (s-i 


) 




beta decay 




15 Mo 


3.39 X 10® 


4.50 X lO'^ 


0.480 


5.17 X 10" 


-7 


"^Fc (29) 


•''•"Fc (25) 


■""'■^Mn (11) 


6.08 X 10" 


11 


^■'Mn (67) 


^^Mn (8) 


32p 


3.82 X 10® 


7.26 X 10^ 


0.464 


3.30 X 10" 


-7 


""Fc (41) 


•"Co (10) 


■''■^Mn (9) 


6.73 X 10" 


-9 


5«Mn (45) 


""Co (18) 


■■■'Mn (11) 


4.13 X 10® 


2.89 X 10* 


0.450 


6.86 X 10" 


-7 


^^Fe (54) 


^^Ni (21) 


Fc (14) 


4.10 X 10" 


-7 


■'«Mn (36) 


•'"V (12) 


"Mn (10) 


4.41 X 10® 


1.30 X 10® 


0.442 


7.57 X 10" 


-6 


■"Fc (22) 


Cr (14) 


^^Mn (13) 


1.74 X 10" 


-6 


"*Mn (34) 


'•'Co (17) 


"Co (12) 


7.25 X 10® 


9.36 X 10® 


0.432 


9.21 X 10" 


-3 


"■'Ni (14) 


•^'®Fc (7) 


52y (.7^ 


8.45 X 10" 


-fi 


"Co (22) 


'•*Mii (19) 


"V (13) 


25 M,,, 


3.79 X 10® 


2.89 X 10^ 


0.487 


3.18 X 10" 


-6 


"''Fe (23) 


""Co (20) 


■'^Ni (19) 


1.53 X 10" 


11 


■'''Mn (49) 


■'■'Fe (17) 


'"^Co (9) 


4.17 X 10® 


3.71 X 10'^ 


0.476 


4.23 X 10" 


-6 


5^Fe (21) 


^^Co (14) 


^^Fe (11) 


8.12 X 10" 


10 


^■^Mn (37) 


5*Co (30) 


^^Fe (8) 


5.03 X 10® 


1.82 X 10* 


0.456 


3.84 X 10" 


-6 


^«Fe (17) 


^^Fe (13) 


^^Ni (10) 


1.00 X 10" 


-6 


5«Mn (45) 


^^V (13) 


«°Co (10) 


5.57 X 10® 


5.05 X 10* 


0.449 


1.45 X 10" 


-5 


5^Fc (16) 


5^Fe (11) 


^^Cr (9) 


7.61 X 10" 


-6 


5'^Mn (19) 


^*Mn (14) 


'^^Cr (10) 


7.75 X 10® 


2.42 X 10® 


0.445 


1.95 X 10" 


-3 


(32) 


^^Cr (9) 


"Fe (7) 


5.17 X 10" 


-5 


^*Mn (18) 


'^'^Cr (13) 


"Mn (7) 



0.450 
0.445 
0.440 



^» 0.435 



0.430 
0.425 
0.420 



2.0 



1.6 
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FIG. 2. Comparison of the center values of Ye (upper part) 
and the iron core sizes (lower part) for 11-40 Mq stars be- 
tween the WW models and the present ones using the LMP 
weak interaction rates. 



the other hand, mainly involves unstable manganese and 
cobalt isotopes. However, the lifetimes are long enough 
to allow for an experimental determination of the rele- 
vant GT strength distribution once radioactive ion beam 
facilities become operational. We finally mention that 
the backresonances contribute noticeably to the stellar 
beta decay rates for these isotopes making also measure- 
ment of the GT_|_ strength on the daughter nuclei (e.g. 
iron and nickel isotopes) very useful. 
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